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Observation of Bound Surface States in Grain Boundary Junctions of High 

Temperature Superconductors 
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We have performed a detailed study of the tunneling spectra of bicrystal grain boundary junc- 
tions (GBJs) fabricated from the high temperature superconductors (HTS) YBa2Cu307_5 (YBCO), 
BiaSraCaCuzOg+i (BSCCO), Lai.85Sro.i5Cu04 (LSCO) and Ndi.85Ceo.i5Cu04-a (NCCO). In all 
experiments the tunneling direction was along the CuOa planes. With the exception of NCCO, for 
all materials a pronounced zero bias conductance peak (ZBCP) was observed which decreases with 
increasing temperature and disappears at the critical temperature. These results can be explained 
by the presence of a dominating d-wave symmetry of the order parameter resulting in the formation 
of zero energy Andreev bound states at surfaces and interfaces of HTS. The absence of a ZBCP 
for NCCO is consistent with a dominating s-wave symmetry of the pair potential in this material. 
The observed nonlinear shift of spectral weight to finite energies by applying a magnetic field is in 
qualitative agreement with recent theoretical predictions. 



There is strong evidence that the superconducting or- 
der parameter (OP) in the HTS has a dominating d-wave 
symmetry [^j^] . For this pairing symmetry there is a tt- 
phase shift of the OP in orthogonal /c-space directions 
resulting in a positive and negative sign of the pair po- 
tential in those directions. This also means that there are 
directions with nodes of the pair potential, e. g. for a pure 
(i2.2_j^2-symmetry, the nodes are along the [110] direction 
in the Cu02 plane. For the tunneling spectra of junctions 
employing HTS electrode materials with a c?-wave sym- 
metry of the OP, a pronounced ZBCP has been predicted 
originating from mid-gap surface (interface) states or zero 
energy bound states (ZES) at the Fermi level The 
physical reason for these states originates from the fact 
that quasiparticles incident and reflecting from the sur- 
face propagate through different order parameter fields 
which leads to Andreev reflection. The constructive in- 
terference between incident and Andreev reflected quasi- 
particles results in bound states. Stable ZES arc formed 
if the scattering induces a change in sign of the OP. For a 
dx2-y2-wa,vc symmetry such sign change and, hence, the 
presence of ZES is possible for all surfaces parallel to the 
c-axis except for those with the lobe directions perpen- 
dicular to the surface, whereas for a s-wave symmetry 
no ZES are possible. The spectral weight of the ZES for 
a (ia;2_j^2-wave symmetry depends on the orientation of 
the surface with respect to the crystal axis. The maxi- 
mum spectral weigth is expected for a (110) surface and, 
hence, a maximum ZBCP is expected for tunneling in 
the direction of the nodal lines, i. e., the [110] direction. 
This has been observed recently using low temperature 
scanning tunneling spectroscopy (LTSTS) j|] and planar 
type junctions [Q. We note that the ZBCP is sensitive 
to surface roughness making it difficult to distinguish be- 
tween the directions in the plane []ll|-[l3|[. 

Initially, the ZBCP in the tunneling spectra of HTS 




junctions has been explained within the Appelbaum- 
Anderson (AA) model |l^ due to the presence of a large 
density of magnetic scattering centers at the surface of 
the junction electrodes. However, the AA-model predicts 
a ZBCP that is not expected to disappear at a certain 
temperature and to split linearly with increasing applied 
magnetic field. Furthermore, it has been suggested re- 
cently that the surface of d-wave superconductors might 
show spontaneously generated surface currents 
and a phase with broken time-reversal symmetry 
In such state the Andreev bound states shift to finite en- 
ergies resulting in a splitting of the ZBCP even in zero 
magnetic field |16|. Applying a magnetic field results in 
a further splitting of the ZBCP. However, in contrast to 
the AA-model prediction this splitting is predicted to in- 
crease nonlinear ly with applied field. In order to clarify 
these issues experimentally and to rule out competing 
explanations for the origin of the ZBCP in GBJs, in this 
report we present a comprehensive analysis of the ZBCP 
for different materials including YBCO (60 K-phase and 
90 K-phase), BSCCO, LSCO, and NCCO. We emphasize 
that for YBCO, BSCCO and LSCO the OP is considered 
to have a dominating c?-wave component [ p^ , whereas for 
NCCO the dominating component most likely is a s-wave 
as suggested by several experiments p^-|2l|]. Therefore, 
if ZES are the origin of the observed ZBCP, such peak 
should be present only for the c?-wave but not for the 
s-wave material. As shown below, for NCCO indeed no 
ZBCP is observed giving strong evidence for the ZES 
scenario and ruling out the magnetic interface scattering 
model. Our data also show a nonlinear evolution of the 
shift of spectral weight to higher energies with increasing 
magnetic field. 

Up to now in several experiments a ZBCP has been 
observed in junctions where only a single electrode was 
based on a cuprate superconductor [p|, p^|j2l| -|25| . Many 
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more experiments with the tunnehng direction along the 
c-axis have been performed, where a ZBCP due to ZES is 
expected only as an artefact of the finite surface rough- 
ness of the HTS electrode. In our experiments we used 
weh defined [001] tilt HTS-GBJs fabricated on bicrys- 
tal substrates It has been shown recently that the 
quasiparticle transport mechanism in these junctions is 
dominated by elastic, resonant tunneling via localized 
states making them suitable for spectroscopic studies 
|p7|-^. A pronounced ZBCP has been observed in the 
tunneling conductance of these GB Js which has been dis- 
cussed both in terms of ZES and the presence of 
magnetic scattering centers at the grain boundary [p9| . 
In this report, we clearly show that the former analy- 
sis can be applied for the HTS-GBJs. There are sev- 
eral advantages of using GBJs. Firstly, these junctions 
are formed by two HTS electrodes and can be fabricated 
easily from different HTS materials [ p6pc| ]. Employing 
mtrinsic interfaces less problems arise from contamina- 
tion due to ex-situ processing of the samples. Secondly, 
the tunneling direction for [001] tilt GBJs is along the 
a6-plane. Thirdly, the direction of tunneling within the 
a6-plane can be varied by varying the misorientation an- 
gle of the bicrystal substrate, although the faceting of 
the grain boundary always results in an averaging over a 
finite range of angles [Q. In this context, we note that 
an exact quantitative description of effects related to the 
faceting is not yet available. 

The GBJs studied in our experiments were prepared 
on symmetrical [001] tilt SrTiOa bicrystals with 24° or 
36.8° misorientation angles. The fabrication and charac- 
terization of the GBJs has been described in detail else- 
where p6| , ^ . The measurements of the current-voltage 
(/(y)) and conductance vs. voltage (G(y) = dI{V)/dV) 
characteristics were performed in a standard four-lead 
arrangement. 
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Differential condy§^^g^ig3(rifi^) voltage for a 36.8° 
[001] tilt LSCO-GBJ between 4 and 23 K. 

Fig. |l| shows a set of typical G(t^)-curves obtained 
for LSCO-GBJs. The critical temperature Tc of the 
LSCO electrodes was about 24 K. Very similar curves (see 
Fig. ^) have been measured for oxygen deficient YBCO 
(Tc « 60 K), BSCCO {Tc « 80 K) and fully oxy- 



genated YBCO (Tc « 90 K). At voltages above the gap 
voltage of the electrode material the G(y)-curves show a 
temperature independent conductance that has an about 
parabolic shape and can be modelled by the influence of 
the applied voltage on the shape of the tunneling barrier 
[ pOf . Below the gap voltage a reduced conductance due 
to a reduced density of states is observed. With increas- 
ing temperature the conductance increases approaching 
the normal state curve GniV) with T approaching Tc. 
The height of the ZBCP decreases with increasing tem- 
perature and vanishes at T = Tc. For most samples 
over a considerable temperature range the decrease fol- 
lows a 1/T dependence. The temperature evolution of 
the G(F)-curves clearly demonstrates that the supercon- 
ducting state is being probed. This proof is important 
with respect to the interpretation of the ZBCP in terms 
of Andreev bound states. The G{V) curves also show 
that the parabohc background conductance G„(V) is a 
normal state effect, which in the following is eliminated 
by normalizing G{V) to G„(F). 




FIG. 2. Normalized conductance vs. normalized voltage 
of [001] tilt GBJs formed by YBCO (90 and 60 K phase), 
BSCCO, and LSCO at T = 4.2 K. In (b) the same depen- 
dence is shown for a NCCO-GBJ. 

In Fig. Ijthe normalized tunneling conductance, G/G„, 
of GBJs fabricated from YBCO, BSCCO, LSCO, and 
NCCO (Tc « 24 K) is plotted versus the voltage normal- 
ized to the gap voltage Vg. Here, Vg — 25, 20, 15, 6, 
and 6 meV was used for BSCCO, YBCO (90 K phase), 
YBCO (60 K phase), LSCO, and NCCO, respectively. In 
first approximation eVg can be considered to be close to 
the gap energy Aq, in contrast to the BCS-theory. See 
for example the calculations in reference ]^ . It is evident 
from Fig. |a that YBCO, BSCCO and LSCO, for which 
the OP is considered to have a dominating d-wave com- 
ponent, qualitatively show the same behavior. A clear 
gap structure with reduced density of states is observed 
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in combination with a ZBCP. For BSCCO and YBCO-90 
the ZBCP is reduced in height as compared to YBCO- 
60 and LSCO. The reason for this reduction is not clear 
at present. However, considering the dependence of the 
ZBCP on the degree of faceting of the grain boundary, 
which determines the amount of averaging over the in- 
plane crystal directions, this observation is not surpris- 
ing. In contrast, for NCCO, which is considered to be a 
s-wave superconductor, only a gap structure but never a 
ZBCP is observed as demonstrated by Fig. ||d. 



in agreement with results published recently pO[ |. 
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FIG. 3. (a) Magnetic field dependence of the ZBCP at 100 
mK for a [001] tilt LSCO-GBJ. In (b), G{V,H) - G{V,0) 
is plotted for the same sample. The applied magnetic fields 
ranged between 1 and 9 T (1.5 T steps). 

We also measured the dependence of the ZBCP on a 
magnetic field H applied parallel to the grain boundary 
plane. A typical result is shown in Fig. ^. The applied 
magnetic field reduces the spectral weight at zero energy 
and shifts it to finite energies that increase with increas- 
ing field. As shown in Fig. this results in a splitted 
peak structure of the difference curve G{V, H) — G{V, 0) 
where the distance between the peaks is defined as 25. 
We emphasize that so far we could not directly observe 
the splitting in the G'(F)-curve down to T = 100 mK, 
where the thermal smearing amounts to only a few 10 
^eV. In Fig. |, 5 is plotted versus H for a LSCO-GBJ 
at 100 mK and a YBCO-GBJ at 4.2 K together with 
data of the direct split in G{H) taken from literature. 
Clearly, 5 does not vary linearly with H as predicted by 
the AA-model. For all investigated samples 5 increases 
slower than linearly and tends to saturate at high fields 
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FIG. 4. Splitting S of thef9(f}^f-^6Ji', 0) curves vs. ap- 
plied magnetic field of a 36.8° [001] tilt LSCO-GBJ at T = 100 
mK and a 24° [001] tih YBCO-GBJ at T = 4.2 K. Also shown 
are data for the direct splitting in G{V, H) of SIN-junctions 
from literature. 

We first will discuss our experimental findings in terms 
of the AA-model jl4) . For tunneling across a barrier con- 
taining localized spins beyond a contribution Gi due to 
direct tunneling without interaction with the spins there 
are two further contributions G2 and G3 to the total con- 
ductivity. The first (G2) is related to tunneling involv- 
ing a spin exchange and the latter (G3) to Kondo-type 
scattering. According to the AA-model |lj] one expects 
G3{V,T) oc ln[EQ/{\eV\ + keT)], where Eg is a cut-off 
energy. Applying a magnetic field the Zeemann splitting 
of the impurity levels causes a dip of width 25 — 2gfiBH 
due to a reduction of G2. Here, g is the factor and /is 
the Bohr magneton. Furthermore, in an applied field the 
Kondo peak is split into three peaks separated by 5 with 
the zero bias peak completely suppressed. This results 
in a peak of G{V, H) - G{V, 0) at eV = ±5. The AA- 
model predicts a peak-to-peak width 25 = 2gfiBH that 
increases oc H, i.e. 5/H = gfiB — 0.12 meV/T for g ~2. 
This is in clear contradiction to our results, which show 
both much larger absolute values of 6/ H up to more than 
2meV/T and a strong increase oi 5/H with decreasing 
H, in agreement with other data reported in literature 
0,11 . Furthermore, for all YBCO, BSCCO and LSCO 
samples the ZBCP always disappeared just at Tc, which 
is significantly different for the different materials. This 
is very difhcult to be explained within the AA-model, 
which predicts the ZBCP to decrease with increasing T 
but not to disappear at a specific temperature. Finally, 
within the AA-model the absence of the ZBCP for NCCO 
would imply the absence of magnetic scatterers for this 
material. Supposing that magnetic scatterers at grain 
boundaries result from oxygen loss and the formation of 
magnetic Cu'^^-ions, the basic difference between NCCO 
and the other materials is difficult to understand. 

We now turn to the Andreev bound state model. As 
discussed above, in this model the ZBCP arises from 
bound states formed by the constructive interference of 
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quasiparticles that propagate through different order pa- 
rameter fieids incident and reffecting from the surface of 
the junction electrode. ZES are formed if the scattering 
induces a change in sign of the OP. Hence, ZES are not 
possible for a s-wave symmetry of the OP. However, in 
the case of a dominating d2.2_j,2 -symmetry of the OP, at 
all surfaces parallel to the c-axis ZES are formed except 
for the surfaces exactly perpendicular to the a- or &-axis 
direction. Hence, the ZES-model naturally accounts for 
the observation that a ZBCP is observed only for YBCO, 
BSCCO and LSCO, which most likely have a dominat- 
ing d-wave component of the OP, whereas it is absent 
for NCCO, which is supposed to have a dominating s- 
wave OP. The ZES-model also qualitatively accounts for 
the increase of the height of the ZBCP with decreasing 
temperature and the nonlinear shift of the peak spectral 
weight to finite voltages with increasing magnetic field 
p3|,pj 33[3^. For example, for a surface to a-axis orien- 
tation of 20°, G{0,T)/Gn{0) was predicted to decrease 
about oc 1/T 1 13 1 in fair agreement with our data. A 
detailed quantitative analysis of our experimental data 
still is not possible, since no prediction of the exact T 
and H dependence of the ZBCP is available taking into 
account the angle averaging due to the faceting of the 
grain boundaries. 

We finally would like to address the possibility of a 
surface state with broken time reversal symmetry as pre- 
dicted by Fogelstrom et al. JM| and experimentally ob- 
served by Covington et al. [p!0|. In this case the ZBCP 
is expected to split in zero magnetic field. Such split- 
ting has not been observed directly in our experiments 
down to temperatures of 100 mK for LSCO and 4.2 K for 
YBCO similar to other experiments |Q. A possible rea- 
son for this observation may be the considerable faceting 
of the grain boundary plane together with impurity scat- 
tering that suppresses the field sphtting of the ZBCP @. 
This is the reason why the observed behavior of S vs. H 
does not provide definitive evidence for a subdominant 
s-wave OP and time-reversal symmetry breaking at the 
grain boundary interface. This issue has to be clarified 
by future experimental and theoretical work taking into 
account the grain boundary faceting and impurity scat- 
tering. 

In conclusion, it has been shown that quasiparticle tun- 
neling in GBJs can be used for probing the symmetry of 
the order parameter in HTS. The tunneling spectra of 
[001] tilt GBJs formed by YBCO, BSCCO and LSCO 
were found to always show a ZBCP while such peak is 
absent for NCCO. The height of the ZBCP decreases 
with increasing temperature and disappears at Tc- These 
observations are not compatible with the assumption of 
tunneling involving magnetic impurities as described by 
the AA-model, but can naturally be explained by the 
presence of zero energy Andreev bound states at surfaces 
of HTS. The existence of ZES represents a further proof 



that the order parameter of YBCO, BSCCO and LSCO 
changes sign on the Fermi surface and most likely has 
a dominating d-wave component. The tunneling data of 
NCCO is consistent with an anisotropic s-wave symmetry 
of the pair potential in the electron doped HTS. The evo- 
lution of the ZBCP with varying applied magnetic field 
and temperature can be qualitatively described within 
the d-wave scenario. 
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